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HIGH  ALTITUDE  ROCKET  PLUME  STRUCTURE: 


EXPERIMENT  AND  CALCULATIONS 


Frederick  P.  Boynton 
Physical  Dynamics,  Inc. 
Detroit,  Michigan 


This  note  presents  some  comparisons  of  local  values 

of  species  number  density  and  impact  pressures  measured  in 

a  simulated  high-altitude  rocket  plume^  with  values  calcu- 

2 

lated  with  a  computer  code  which  the  author  has  used  ex¬ 
tensively  over  the  past  several  years.  The  measurements 
were  conducted  at  the  Von  Karman  Gas  Dynamics  Facility 

( VKF )  of  the  U.S.  Air  Force  Arnold  Engineering  Development 

.  2 
Center  and  included  pitot-tube  determinations  of  pu  and 

electron-beam  fluorescence  measurements  of  the  local  num¬ 
ber  density  n^  of  plume  and  free-stream  gases. 

Interpretation  of  the  measurements  at  the  experimental 
conditions  is  not  straightforward,  since  significant  correc¬ 
tions  to  the  raw  data  are  required  in  order  to  determine  ab¬ 
solute  values  of  pu2  and  n..  The  author  has  discussed  the 
experimental  results  with  several  people  connected  with 
the  experiment.*  Their  consensus  is  that,  the  absolute  val¬ 
ues  of  the  data  of  Reference  )  may  require  further  correction. 

*  The  author  particularly  wishes  to  acknowledge  discussions 
with  Drs.  L.  Ouinn,  Air  Force  Rocket  Propulsion  Laboratory; 
J.D.  Stewart,  Aerospace  Covp.j  and  E.  A.  Sutton,  Aerodyne 
Research,  Inc. 
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Thus,  at  this  time  only  a  semi-quantitative  comparison  with 
the  code's  predictions  is  possible.  On  the  other  hand,  these 
data  represent  the  only  available  experiments  on  local  plume 
conditions  in  the  range  of  shock  layer  Reynolds  number  char¬ 
acteristic  of  large  vehicle  plumes  between  100  and  200  km 
altitude,  which  is  the  regime  for  which  the  code  was  origi¬ 
nally  intended.  A  comparison  between  the  data  and  a  calcu¬ 
lation  should  be  useful,  even  though  some  further  data  re¬ 
duction  may  be  needed  to  make  the  comparison  completely  valid. 

For  the  comparison,  we  choose  test  condition  IC5,  con¬ 
ducted  in  Tunnel  M  at  VKF .  In  this  test  a  helium  plume  was 
released  from  a  5  psia  chamber  through  a  conical  nozzle  of 
1.61  area  ratio  into  a  Mach  18.15  nitrogen-free  stream  at 
6y  Hg  static  pressure.  On  Simons'  map  of  plume  regimes, 
this  flow  field  would  lie  slightly  above  the  lines  correspond¬ 
ing  to  merging  of  the  barrel  shock  and  the  low-altitude  boun¬ 
dary  of  the  transitional  regime.  The  computer  code  performs 
a  "viscous  layer"  calculation  in  which  the  shocks  are  assumed 
to  be  thin.  Between  the  two  shocks,  the  equations  of  motion 
consist  of  the  full  Euler  equations  with  viscous  terms  resem¬ 
bling  those  in  the  boundary  layer  equations.  (See  Reference  2 
for  a  complete  description.)  The  mixing  zone  between  plume  and 
free-stream  gas  need  not  be  thin  with  respect  to  the  shock 
layer  as  a  whole.  We  should  therefore  expect  that  this  case 


*  A  comprehensive  review  of  the  data  is  currently  being  con¬ 
ducted  by  Aerodyne  Research,  Inc.,  and  it  is  hoped  that 
either  reliable  values  or  reliable  error  estimates  will  be 
available  within  the  next  year. 
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should  provide  a  stringent  test  of  the  calculation's  ability 
to  describe  conditions  in  the  mixing  region  between  the 
plume  and  free  stream. 


conditions 

;  of  the  experiment  and  calculation 

are 

Table  I. 

For  the  calculation 

it  was  assumed 

that 

/ 

TABLE  I 

Exhaust 

and  Free  Stream  Tost 

Conditions 

Exhaust  Gas 

Free  Stream 

P 

o 

258  mm  H 

g 

To 

7  70  0  K 

2  900  °  K 

p 

28.1  mm  Hg 

6.00  Hg 

T 

31 7  °K 

X 

O 

O 

in 

u 

2.17  *  105  cm/sec 

2.49  *  105  cm/sec 

M 

2.07 

18.15 

gas 

He 

N2 

Y 

1.667 

1.4 
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all  mixing  was  steady  and  laminar,  that  the  viscosities  of 
N2  and  lie  could  be  adequately  represented  as  power-law  func¬ 
tions  of  temperature,  and  that  Prandtl  and  Schmidt  numbers 
wore  constant.  The  experimental  setup  includes  a  forebode 
needed  to  house  the  nozzle;  this  forobody  was  a  core-cylinder 
with  a  chamfered  boat-tail.  The  inviscid  flow  over  this  fore¬ 
body  was  calculated  and  included  in  the  description  of  the 
free  stream.  Any  viscous  effects  due  to  boundary  layer  for¬ 
mation  or  separation  are  not  included  in  the  calculation. 
Except  for  the  nose  region  of  the  plume,  where  the  air  shock 
is'  detached,  the  calculations  were  performed  with  the  MUL7I- 

TUBi:  code.  In  the  nose  region,  a  recently  developed  thin 
4 

layer  code  was  employed.  Partly  because  of  the  configura¬ 
tion  of  the  nose  region  of  this  plume,  this  calculation's 
description  of  the  nose  region  is  poor,  and  it  should  be 
regarded  more  as  a  means  of  providing  reasonable  initial 
conditions  to  the  wholly  supersonic  part  of  the  plume  than 
as  a  realistic  description  of  the  nose  region.  That  the 
initial  conditions  are  in  fact  adequate  is  evidenced  by  the 
behavior  of  the  supersonic  calculations,  which  exhibit  only 
mild  pressure  oscillations  in  the  first  centimeter  or  two 
downstream  of  the  transition  between  the  two  methods  of  cal¬ 
culation.  The  results  of  the  calculation  in  the  forward  re¬ 
gion  of  the  plume,  including  the  forebody  flow,  are  shown  in 
Figure  2. 

Impact  pressures  as  a  function  of  axial  distance  from 
the  nozzle  are  shown  in  Figure  3.  The  measured  values  lie 
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close  to,  but  slightly  below,  the  computed  values  until  about 
15  cm  downstream,  at  which  point  the  measured  values  rise 
above  the  computed  values  and  actually  increase  beyond  35  cm 
downstream.  This  nest  likely  reflects  the  upstream  influ¬ 
ence  of  the  Mach  disc,  which  is  expected  to  be  very  diffuse 
in  this  plume. 

Inpact  pressures  as  a  function  of  radial  distance  from 
the  centerline  are  shown  in  Figures  4-6  at  distances  of  3, 

5,  and  10  inches  (7.6,  12.7,  and  25.4  cm)  downstream.  The 
major  apparent  qualitative  difference  between  calculation 
and  experiment  is  the  thickness  of  the  experimental  shocks 
(expected  on  the  basis  cf  Simons'  regime  map).  Absolute 
values  of.  the  measured  and  calculated  ou  agree  rather  well 
in  the  internal  region  of  the  shock  layer  and  somewhat  less 
well  elsewhere.  (Since  further  correction  to  the  measurements 
ma /  be  needed,  one  cannot  make  too  much  of  any  absolute  com¬ 
parisons  at  present.)  It  docs  appear  from  the  measurements 
that  the  jot  shock  is  quite  thick.  The  predicted  peak  in 
the  impact  pressure  profiles  on  the  inner  side  of  the  plume 
shock  layer  (a  result  of  competition  between  an  outward  pres¬ 
sure  gradient  duo  to  centrifugal  forces  and  an  outward  tem¬ 
perature  gradient  due  to  heat  conduction  from  the  shocked  N2) 
is  only  found  at  the  three-inch  station.  Outside  of  these 
thick  shocks,  the  agreement  of  calculation  and  experiment  is 
quite  satisfactory. 
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Mo  and  molecule  number  densities  determined  fro n 

el octron-boan  fluorescence  measurements  are  compared  wi*-h 

the  calculated  values  in  Figures  7  -  10.  One  should  be 

aware  that  at  the  time  these  measurements  were  made  the 

electron-bean  technique  va.i»  still  under  development  at  VKF,  and 

that  many  questions  remain  concerning  absolute  calibrations 

in  llo-N^  mixtures,  including  beam  spread,  secondary  electron 

effects,  and  quenching  of  excited  states.  Where  the  flow  is 

mostly  N2,  a  good  check  is  afforded  by  comparison  with  the 

known  f roc-stream  conditions.  The  values  of  n,,  quoted  in 

no 

Reference  1  are  consistently  40  -  60*  larger  than  those  de¬ 
termined  from  the  calculation.  If  one  accepts  the  quoted 
nozzle  properties  (from  which  one  can  calculate  the  num¬ 
ber  flow  of  Me  atoms),  then  cither  the  calculated  velocities 
in  the  Me  plume  are  very  badly  wrong,  a  circumstance  which 
would  do  violence  to  a  number  of  long-accepted  concepts  in 
gas  dynamics,  or  there  is  a  sizeable  correction  required  in 
the  reported  absolute  He  number  densities. 

The  qualitative  behavior  of  the  calculated  and  measured 
values  of  is  substantially  the  same  except  in  the  region 

of  the  jot  shock.  In  fact,  simply  reducing  the  reported 
values  of  nHe  by  40%  brings  them  into  rather  satisfactory 
quantitative  agreement  with  the  calculation.  Conservation  of 
helium  molecules  at  velocities  approaching  the  limiting  velo¬ 
city  (2.5  to  2.8  km/sec)  requires  a  correction  factor  of  this 
order,  though  of  course  it  may  be  different  in  different  re- 
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qions  depending  upon  density,  temperature,  and  composition. 


The  calculated  and  measured  values  of  n^^.'  excepting 
the  region  of  the  free-stream  shock,  are  in  very  satisfac¬ 
tory  agreement  at  the  10-inch  station  (Figure  9).  However, 
at  locations  further  upstream  the  measurements  show  consid¬ 
erably  greater  penetration  of  ^  molecules  at  the  inner 
edge  of  che  jet  shock.  This  anomalously  high  diffusion  is 

also  evident  in  the  results  of  test  IC2  of  this  series  and 

5 

in  several  tests  conducted  in  the  10-V  tunnel.  The  dis¬ 
crepancy  between  the  observed  and  predicted  diffusion  of 
into  the  jet  is  disturbing,  since  it  suggests  that  tempera¬ 
ture  profiles  could  also  be  broader  than  calculated.  If 
this  is  the  case,  estimates  of  emission  from  the  forward 
part  of  the  plume  may  be  significantly  in  error. 

Several  possible  reasons  for  the  discrepancy  have  deen 
examined,  and  none  seems  to  supply  a  satisfactory  explanation 
of  it.  It  is  difficult  to  ascribe  it  to  some  rarefied  flow 
phenomeron,  since  an  average  Nj  molecule  entering  the  plume 
upstream  from  inis  region  suffers  10  to  20  collisions  before 
getting  across  the  shock  layer.  Thus  the  proper  description 
of  molecular  motion  is  in  terms  of  diffusion  rather  than 
molecular  penetration,  and  this  description  is  embodied  in 
our  calculation.  It  is  possible  that  the  forebody  introduces 
some  additional  mixing  into  the  forward  part  of  the  plume. 

We  cannot  assess  this  effect  with  available  codes.  Discus- 
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sions  with  people  responsible  for  the  tests  indicate  that 
great  care  was  taken  to  align  the  electron  beam  in  the 
Tunnel  M  test  series,  and  it  is  doubtful  that  measurement 
errors  arise  from  that  source.  There  are  other  possible 
errors,  such  as  those  due  to  beam  spread,  secondary  electron 
effects,  and  movement  of  the  molecules  before  emission. 
Assessment  of  these  effects  is  currently  under  way  else¬ 
where,  and  we  should  properly  reserve  comment  until  this  re¬ 
interpretation  of  the  data  is  complete.  However,  we 
note  that  there  is  a  difference  in  radial  location  of  the 
measured  maxima  in  pu2  and  n^  for  the  3-inch  and  5-inch  sta¬ 
tions,  while  we  should  expect  these  to  be  coincident. 

In  summary,  the  available  data  provide  only  a  semi- 
quantitative  test  of  the  computational  procedure  because  of 
questions  regarding  calibration.  Within  these  limits,  the 
agreement  of  experiment  and  calculation  within  most  of  the 
mixing  region  between  plume  and  free  stream  is  acceptable. 

The  apparent  deep  penetration  of  N2  molecules  into  the  plume 
near  the  nose  is  unexplained.  Unfortunately,  no  comparison 
can  currently  be  made  with  local  temperatures,  which  are  of 
considerable  importance  to  radiation  calculations. 
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Figure  1.  Sketch  of  the  Flow  Field  of  a  Rocket  Exhaust 
Plume  at  High  Altitudes  (continuum  flow) 


Figure  2.  Calculated  Structure  of  the  Flow  for  VKF  Test  Condition  IC5, 
Tunnel  V\,  in  the  Forward  Portion  of  the  Plume.  Flow  over  body 
CALCULATED  WITH  MULTITUDE;  FLOW  IN  THE  PLUME  SHOCK  LAYER  TO 
JUNCTION  LINE  CALCULATED  WITH  THIN  LAYER  CODE;  REMAINDER  OF  FLOW 
CALCULATED  WITH  MULTITUDE.  SHOCK  LAYER  IS  VISCOUS;  AND  DIVIDING 
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Figure  3.  Impact  Pressure  Along  Plume  Centerline  as  a 
Function  or  Distance  from  Nozzle 
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figure  4.  Impact  Pressure  vs  Radial  Distance  in  a  Plane 
3  Inches  Downstream  from  Nozzle  Exit 
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t  igure  5.  Impact  Pressure  vs.  Radial  Di stance  in  a  Plane 
5  Inches  Uownstream  from  Nozzle  Exit 
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Figure  6.  Impact  Pressure  vs.  Radial  Distance  in  a 
Plane  10  Inches  from  Nozzle  Exit 
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Figure  7,  Number  Densities  of  He  and  N2  vs.  Radial 
Distance  in  a  Plane  3  Inches  Downstream 
from  Nozzle  Exit 
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Figure  3.  Number  Densities  of  He  and  N2  vs.  Radial 
Distance  in  a  Plane  5  Inches  Downstream 
from  Nozzle  Exit 
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Figure  9,  Number  Densities  of  He  and  N2  vs.  Radial 
Distance  in  a  Plane  10  Inches  Downstream 
from  Nozzle  Exit 
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